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ABSTRACT 

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. 
We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - transition of OH + and the 1115 GHz 
In ~ Ooo transition of ortho-H 2 + . The resultant spectra show absorption by ortho-H 2 + , and strong absorption by OH + , in foreground material 
at velocities in the range to 70km s with respect to the local standard of rest. The inferred OH + /H20 + abundance ratio ranges from ~ 3 to 
~ 15, implying that the observed OH + arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the 
distribution of OH + and H20 + in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption 
measurements, and dissimilar from that typical of other molecular tracers. The observed OH + /H abundance ratio of a few xl0~ 8 suggests acosmic 
ray ionization rate for atomic hydrogen of 0.6 - 2.4 x 10~' 6 s~', in good agreement with estimates inferred previously for diffuse clouds in the 
Galactic disk from observations of interstellar Hj and other species. 

Key words. ISM: molecules - Submillimeter: ISM - Astrochemistry - Molecular processes 
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1. Introduction 

After hydrogen and helium, oxygen is the most abundant heavy 
element in the Universe, accounting for more than 1% of the 
mass of baryons in the Galaxy. Of the ~ 150 distinct molecules 
detected in the interstellar gas to date, more than one-quarter 
contain oxygen, including many of the most widely-observed 
species such as CO, OH, H z O, CH 3 OH, H 2 CO, SiO, and S0 2 . 
Not surprisingly, the chemical processes leading to the in- 
corporation of interstellar oxygen atoms into oxygen-bearing 
molecules have been the subject of numerous theoretical inves- 
tigations. These studies have identified three types of process 
that can be important in various interstellar environments: (1) in 
warm regions, where the gas temperature exceeds ~ 300 K as 
a result of heating by shocks, by turbulent dissipation, or by in- 
frared or ultraviolet radiation, the oxygen chemistry can be initi- 
ated by the endothermic reaction of atomic oxygen with H 2 ; (2) 
on grain surfaces, the formation of oxygen-bearing molecules 
can occur following the adsorption of atomic oxygen; the resul- 
tant molecules can then be released into the gas-phase by pho- 
todesorption or evaporation if the grains are exposed to ultravi- 
olet radiation or heated sufficiently; and (3) in gas clouds that 
are irradiated by cosmic rays or X-rays, the oxygen chemistry is 
controlled by a series of ion-neutral reactions, and is initiated by 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



the formation of the hydroxyl cation OH + , which can react with 
molecular hydrogen to form H 2 + and then H 3 + . 

This third type of process, ion-neutral chemistry, is believed 
to dominate the formation of oxygen-bearing molecules within 
the cold quiescent interstellar medium (ISM). However, the de- 
tection of two key intermediaries in this reaction sequence, the 
OH + and H 2 + molecular ions, has proven elusive. Like other 
light hydrides with small moments of inertia, these molecules 
have rotational transitions at high frequencies that are difficult 
or impossible to observe from ground-based observatories. Until 
very recently, neither OH + nor H 2 + had been detected in the 
interstellar gas (although optical transitions of H 2 + have long 
been detected in comets, e.g. Wehinger et al. 1974). In the last 
three months, however, the observational picture has changed 
rapidly. With the availability of the HIFI instrument (de Graauw 
et al. 2010) on Herschel (Pilbratt et al. 2010), detections of 
H 2 + have been reported in absorption toward DR21, Sgr B2 
(M), NGC6334I (Ossenkopf et al. 2010) and G10.6 - 0.4 (Gerin 
et al. 2010). In G10. 6-0.4, which was targeted for absorption 
line studies as part of the PRISMAS ("PRobing InterStellar 
Molecules with Absorption line Studies") key program, absorp- 
tion by H 3 + and strong absorption by OH + were also reported. 
In addition, Her.sc/ie//SPIRE observations of Mrk 23 1 have re- 
vealed luminous OH + emission (van der Werf et al. 2010), while 
recent ground-based observations have led to detection of inter- 
stellar OH + in the near-UV spectral region (Krelowski, Beletsky 
& Galazutdinov 2010; who used the ESO Paranal observatory 
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Fig. 1. Spectra of H z O + In - o (green) and OH + N = 1 - 
(blue) transitions obtained toward W49N. The velocity scale ap- 
plies to the strongest hyperfine component, with assumed fre- 
quencies of 971803.8 MHz for OH + and 1115204 MHz for 
H20 + . Note that because HIFI employs double sideband re- 
ceivers, the complete absorption of radiation at a single fre- 
quency will reduce the measured antenna temperature to one- 
half the apparent continuum level. 



to detect a weak absorption line at 358.3769 nm) and at submil- 
limeter wavelengths (Wyrowski et al. 2010, who used the APEX 
observatory to detect OH + in absorption toward Sgr B2). 

In this Letter, we report Herschel/HIFI observations of OH + 
and H20 + towards a second strong continuum source that we 
have targeted in the PRISMAS program, W49N. This luminous 
region of star formation, lying within the W49A complex, is lo- 
cated at a distance of 11.4 kpc, a value reliably determined by 
means of maser proper motion studies (Gwinn et al. 1992). The 
sight-line to W49N is known to intersect a large collection of 
foreground gas clouds, with velocities relative to the local stan- 
dard of rest (LSR) ranging from ~ to ~ 70 km s , these have 
been widely detected by means of absorption line spectroscopy 
of multiple species, including H (Brogan & Troland 2001, Fish 
et al. 2003), H 2 0, OH, (Plume et al. 2004), O (Vastel et al. 2000), 
HCO\ HCN, HNC, and CN (Godard et al. 2010). 

2. Observations and data reduction 

In the observations reported here, we targeted two transitions - 
the 972 GHz N = 1 - transition of OH\ and the 1 1 15 GHz 
In - Ooo transition of ortho-H20 + - in the lower sidebands of 
Bands 4a and 5a of the HIFI receiver. To help determine whether 
any observed feature did indeed lie in the expected sideband, 
three observations were carried out with slightly different set- 
tings of the local oscillator (LO) frequency. These observations, 
with on-source integration times of 59 s (OH + line) and 89 s 
(H20 + line) each, were carried out on 2010 April 18, using the 
dual beam switch (DBS) mode and the Wide Band Spectrometer 
(WBS). The WBS has a spectral resolution of 1.1 MHz, corre- 
sponding to a velocity resolutions of 0.34 km s -1 to 0.30 km s~' 
at the frequencies of the OH + and H20 + transitions. The tele- 
scope beam was centered at a =19hl0ml3.2s, 6 - 09°06'12.0" 
(J2000). The reference positions for these observations were lo- 
cated 3' on either side of the source along an East- West axis. 

The spectroscopic parameters for the observed transitions 
have been summarized by Ossenkopf et al. (2010) and Gerin 
et al. (2010). Both the OH + and H 2 + transitions exhibit hy- 
perfine structure associated with the interaction of the nuclear 



spin magnetic moment with that resulting from electronic spin. 
The OH + N - 1 —0 transition has three hyperfine components, 
at velocities of -35.6, -0.5, and km s relative to the cen- 
tral hyperfine component, and with absorption line strengths in 
the ratio 1 : 5 : 9 if the lower levels are populated in propor- 
tion to their statistical weights. For H20 + , there are five hyper- 
fine components; the corresponding velocity shifts are -25.5, 
-15.9, 0, 4.8 and 14.4 km s and the ratio of line strengths is 
1 : 8 : 27 : 8: 10 . Frequencies accurate to 1.5 MHz or better are 
available for the OH + transition, but those for the H20 + transi- 
tion have been the subject of some controversy (Ossenkopf et al. 
2010). As discussed in Appendix A, we favor a frequency close 
to that given by Miirtz et al. (1998) for the strongest hyperfine 
component, 1115.204 + 0.002 GHz. 

The data were processed using the standard HIFI pipeline 
to Level 2, providing fully calibrated spectra of the source. The 
Level 2 data were analysed further using the Herschel Interactive 
Processing Environment (HIPE; Ott 2010), version 2.4, along 
with ancillary IDL routines that we have developed. For each 
of the target lines, the signals measured in the two orthogonal 
polarizations were in excellent agreement, as were spectra ob- 
tained at the three LO settings when assigned to the expected 
sideband. We combined the data from the three observations, 
and from both polarizations, to obtain an average spectrum. 

3. Results 

Figure 1 shows the WBS spectra for the OH + N - I - 
(blue) and ortho-H20 + In - Ooo (green) transitions, with the 
frequency scale expressed as Doppler velocities relative to the 
Local Standard of Rest (LSR) for the strongest hyperfine com- 
ponent. The double sideband continuum antenna temperatures 
are ^(cont) = 7.89 K and 9.24 K respectively for the OH + and 
ortho-H20 + transition, and the r.m.s noise values are 0.069 K 
and 0.1 17 K. Because HIFI employs double sideband receivers, 
the complete absorption of radiation at a single frequency will 
reduce the measured antenna temperature to rougly one-half the 
apparent continuum level (given a sideband gain ratio ~ unity). 
Horizontal lines in Figure 1 indicate the values of ^(cont) and 
0.5 ^(cont). The OH + absorption feature shows a flat bottom in 
the LSR velocity intervals [35,42] and [60,70] km s -1 that does 
not exactly fall to 0.5 ^(cont). This behavior suggests that the 
sideband gain ratio (SBR) is not exactly unity, and we therefore 
treat the SBR as an adjustable parameter in the line fitting pro- 
cedure described below. 

In fitting the observed spectra, we have assumed the absorp- 
tion to arise in a set of clouds with Gaussian opacity profiles. 
Moreover, given the close chemical relationship between the two 
cations that we have detected, we assumed that the OH + and 
ortho-H20 + arise within the same set of clouds. Each cloud is 
therefore characterized by a set of four quantities that we treat 
as adjustable parameters: the centroid velocity, \lsr, the veloc- 
ity width, Av (defined here as the full-width-at-half-maximum), 
and the OH + and H20 + column densities. Here, we neglected 
any contribution from para-H20 + to the total H20 + column den- 
sity; the loi - loi transition of para-H20 + was searched for but 
not detected toward W49, and observations of para-H20 + toward 
Sgr B2 (M), reported recently by Schilke et al. (2010), imply that 
the ortho-to-para ratio is typically large (~ 5) in diffuse clouds 
along that sight-line. 

We also assumed that the population of OH + and ortho-H20 + 
in excited rotational states is negligible, and that the individual 
hyperfine states within the ground rotational state are populated 
in proportion to their statistical weights. Our line fitting proce- 
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Table 1. W49N absorption line components 



VLSR 


Av 


W(OH + ) 


JV(H 2 + ) 


OH + /H 2 + 


/(H 2 )° 


kms~' 


kms~' 


10 13 cnr 2 


10 13 cnr 2 






6.6 


5.1 


2.9 


0.42 


6.8 


0.034 


13.4 


4.6 


2.2 


0.69 


3.2 


0.082 


21.2 


7.0 


8.2 


0.56 


14.6 


0.015 


34.7 


9.0 


26.1 


2.68 


9.7 


0.023 


43.5 


8.4 


7.5 


0.72 


10.5 


0.021 


62.5 


13.3 


18.9 


1.67 


11.3 


0.019 



"using eqn. (1) for an assumed temperature of 100 K and an assumed 
electron abundance of 1.4 x 10~ 4 . The required molecular fraction, 
/(H 2 ), scales as T~° 5 x e 



dure included two additional free parameters: (1) the SBR for the 
OH + observations; and (2) the rest frequency of the H20 + tran- 
sition (see Appendix A) We adopted the minimum number of 
absorption components needed to provide a reasonable fit to the 
data - six - and obtained initial estimates of the cloud velocities 
and FWHM's by eye. These estimates were then refined by us- 
ing the IDL mpfitfun routine, which implements the Levenberg- 
Marquardt method for obtaining a least-squares fit. 

Figure 2 shows the best fit that could be obtained to the OH + 
(top panel) and H20 + (middle panel) spectra for a model with 
six foreground clouds. The reduced^ 2 for this fit was 1.83, sug- 
gesting that the "minimal" absorption model does not entirely 
account for the data. The best-fit cloud parameters are given 
in Table 1 . The best-fit line widths are moderately broad, par- 
ticularly for the 62.5 km s _1 feature, and might - in reality - 
represent the superposition of multiple narrower features. The 
SBR derived for the OH + spectrum was 0.884 (signal sideband 
gain divided by image sideband gain), and the frequency derived 
for the strongest H20 + hyperfine component was 1 1 15209 MHz 
(see Appendix A). In the bottom panel of Figure 2, we present 
the OH + (blue) and H20 + column densities (green) obtained per 
unit velocity width in our best fit model. The cloud parameters 
in Table 1 appear to be relatively robust. Although the OH + ab- 
sorption is completely thick for LSR velocities in the ranges 30 - 
42 and 60 - 70 km s -1 (as determined for the strongest hyperfine 
component), the existence of a weaker hyperfine component at 
an offset of -35 km s~ ] means that arbitrarily large column den- 
sities of OH + can only be accommodated within a narrow range 
of LSR velocities: 65 - 70 km s" 1 . 

4. Discussion 

For diffuse molecular clouds, theoretical models (e.g. 
van Dishoeck & Black 1986) have elucidated the ion-neutral 
chemistry that leads to the production of oxygen-bearing 
molecules. The reaction network is initiated by the cosmic -ray 
ionization of atomic hydrogen to form H + , which undergoes 
charge transfer with atomic oxygen in a reaction that is slightly 
endothermic by the equivalent of ~ 230 K. The resultant + 
ions can then undergo a series of exothermic hydrogen atom 
abstraction reactions, leading to the OFF, H20 + and H30 + . The 
latter is removed by dissociative recombination, which can be 
the dominant source of the neutral molecules OH and H2O. If 
the ratio of electron density to H2 is sufficiently large, as it is in 
diffuse clouds of small molecular fraction, the pipeline leading 
from + to OH + to H20 + to H30 + can be a leaky one, with 
the flow of ionization reduced at each step by the dissociative 
recombination of OH + and H20 + . In dense molecular clouds, 
by contrast, the atomic hydrogen abundance and temperature 
are both too small for + production to be efficient, and the 
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Fig. 2. Fit to the OH + (top panel) and H20 + (middle panel) ab- 
sorption spectra observed toward W49N. The black curve shows 
the data, and the red curve is the fit obtained for the cloud pa- 
rameters listed in Table 1 . The corresponding column densities 
per unit velocity interval are shown in the bottom panel, for OH + 
(blue) and H20 + (green). Also shown, for comparison, are the HI 
(black; from Fish et al. 2003) and para-H20 (red; Sonnentrucker 
et al. 2010) column densities. Where the HI column density ex- 
ceeds the maximum value on the vertical scale, that value should 
be regarded as a lower limit, the 21 cm transmission being zero 
to within the errors. 



dominant source of OH + is the reaction of H3 with O. Once 
again, the chemistry is driven by cosmic ray ionization - the 
original source of H^ - but now the conversion of OH + to H20 + 
to H30 + proceeds with almost 100% efficiency. 

As discussed by Gerin et al. (2010), the OH + /H 2 + abun- 
dance ratio provides a critical probe of the molecular frac- 
tion. Regardless of the production mechanism for OH + , the 
OH + /H20 + ratio is given by 



ft(H 2 |OH + ) »(e)fc(e|H + ) 

k(H 2 \H 2 + ) + n(H 2 )/t(H 2 |H 2 + ) 



= 0.64 + 1490 ■ 



XpL -, 



/(H 2 ) 



(1) 



where k(X\Y) is the rate coefficient for reaction of X with Y, 
n(X) is the density of species X, x e - «(e)/[2n(H2) + «(H)] is 
the fractional ionization, /(H2) = 2«(H2)/[2«(H2) + «(H)] is the 
molecular fraction, and T = 100 T 2 K is the temperature. Here, 
we adopt the reaction rate coefficients tabulated in Gerin et al. 
(2010; their Table 3). 

This simple expression, which accounts exactly for all re- 
sults that we have obtained to date with the Meudon PDR model 
(Le Petit et al. 2006; Goicoechea & Le Bourlot 2007), assumes 
only that H20 + is produced by reaction of OH + with H2, is 
destroyed by dissociative recombination or reaction with H2, 
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Table 2. OH + and H 2 + abundances 



VLSR N(H) a 

krns" 1 10 21 cnT 2 



OH+/H H 2 Q + /H f H e/«2 



30-50 

50-78 



6.95 

7.23 



5.1 x 10- s 
3.0 x 10" 



5.1 x 10^ 
2.6 x 10~ 9 



1.15 
0.63 



"from Godard et al. (2010) for an assumed spin temperature of 100 K. 

and has reached a steady-state abundance. For the six absorp- 
tion components in which the OH + and H20 + column densities 
are well determined, the resultant abundance ratios range from 
~ 3 — 15, requiring /(H2) in the range ~ 0.02 - 0.08. Since the 
largest plausible fractional ionization in a neutral gas cloud is 
~ 1.4 x 10~ 4 , corresponding to the complete ionization of car- 
bon, the OH + ions must reside primarily within clouds of low 
molecular fraction. The derived molecular fraction is given in 
Table 2 for each absorption component. 

Our conclusion that OH + resides primarily in clouds that are 
predominantly atomic - a result obtained previously by Gerin et 
al. (2010) for the sight-line to G 10.6-0.4 - is strongly corrob- 
orated by the observed distribution of absorbing material in ve- 
locity space (Fig. 2, bottom panel). The distribution of H20 + (in 
green) and OH + (blue) is similar to that of atomic hydrogen, as 
determined by 21 cm observations reported by Fish et al. (2003; 
black), (although atomic hydrogen is evidently present without 
accompanying OH + and H20 + at LSR velocities in the range 
-10 to km s _1 ). The distribution of OH + and H 2 + is strik- 
ingly dissimilar from that of para-F^O (red), which has been 
observed by HIFI and reported by Sonnentrucker et al. (2010). 
Whereas water and other molecules detected previously toward 
W49N appear to reside primarily in clouds of relatively narrow 
velocity width, OH + , Fi20 + and H are more broadly distributed 
in velocity space. 

In Table 2, we present the abundances of OFF and H20 + rel- 
ative to atomic hydrogen in two velocity ranges considered pre- 
viously by Godard et al. (2010). Here we adopt the estimates of 
7V(H) derived by Godard et al. from the HI 21 cm spectra of Fish 
et al. (2003) for an assumed spin temperature, 7$, of 100 K. As 
discussed by Gerin et al. (2010), the OH + to atomic hydrogen ra- 
tio probes the cosmic ray ionization rate, (u, defined here as the 
total rate of ionization per hydrogen atom. The presence of OH + 
in clouds of small molecular fraction was anticipated by Liszt 
(2007), in his investigation of the time-scale for H2 formation 
in diffuse atomic clouds. In his models, the hydrogen molecular 
fraction approaches steady state only after 10 7 yr. His predicted 
formation rate of OH + at £# = 3 x 10~ 16 s" 1 is sufficient to attain 
n(OH + )/n H - 5 x 10~ 8 in clouds at n H ~ 50 cirT 3 when the 
cloud age is in the range 2 to 10 million yr. 

Defining e as the ratio of the OH + production rate to the cos- 
mic ray ionization rate, assuming that OH + is destroyed by re- 
action with H2 and dissociative recombination at a rate equal to 
its formation rate, and using eqn. (1) to relate n e ln(Ylz) to the 
observed OH + /H20 + ratio, we find that 



Cne ___ 17 *(OH + )r-°- 5 



= 9.1x10" 



"2 



10- 



1.410" 



l + ^_] s -i, (2) 

R - 0.64 ' 



where R = «(OH + )/«(H 2 + ), *(OH + ) = n(OH + )/n H , n 2 = 
7?H/10 2 cirT 3 and Hh - «(H) + 2n(H2). Values for £#e/«2 are 
tabulated in Table 2 for an assumed temperature of 100 K and a 
fractional ionization of 1.4 \0 A . If cosmic ray ionization is en- 
tirely responsible for the production of OH + , then e represents 
the efficiency with which ionization is transferred from H + to 



OH + . In this case, e must be < 1, and our analysis yields a ro- 
bust lower limit on £#. 

Based upon parameter studies performed with the Meudon 
PDR model, to be described in detail in a future publication, 
we find that values of e close to unity can be achieved under 
a wide range of conditions typical of diffuse molecular clouds. 
We have computed the structure of diffuse molecular clouds for 
a set of models with gas-phase carbon and oxygen abundances 
of 1.38 x 10~ 4 and 3.02 x 10~ 4 respectively with respect to H nu- 
clei, cosmic ray ionization rates £h of 10~ 17 , 10~ 16 , or 10~ 15 s , 
density «h of 10 2 or 10 3 cirr 3 , and UV field xuv of 1, 10, or 
100 (normalized relative to the mean interstellar value). Under 
most of the conditions that we considered, e lies in the range 0.5 
- 1.0 whenever R lies in the observed range of 3 - 15, the only 
exceptions being cases with (V) Xuv - 100 and £# = 10 -17 s _1 , 
for which e can reach values as large as 3; and (2) xuv - 1 and 
«h = 10 3 cm -3 , for which e lies in the range 0.02 - 0.1. In the 
first of these exceptional cases, the temperature reaches several 
hundred Kelvin and the production of OH + is enhanced by reac- 
tion of O with H2 to form OH, followed by photoionization of 
OH to produce OH + ; in the second case, the temperature is too 
low to permit efficient charge exchange between O and H + . 

Given the temperature, 100 K, and density, 100 cirT 3 , typical 
of diffuse clouds with a small molecular fraction, our results sug- 
gest that the OH + and H20 + abundances observed toward W49 
imply a cosmic ray ionization rate £3 ~ 0.6 - 2.4 x 10 _I6 s _I 
for atomic hydrogen. This is the total rate of ionization per hy- 
drogen atom, including secondary ionizations. The range given 
here reflects uncertainties in e, assumed to lie between 0.5 and 
1 , and variations in R from one component to another; for tem- 
peratures and densities different from those assumed above, the 
derived value of (h would scale as nT~ 05 T s l . This estimate is 
in good agreement with values for the primary ionization rate 
per H nucleon obtained by Indriolo et al. (2007) from an analy- 
sis of HJ abundances observed along many sight-lines intersect- 
ing diffuse clouds in the Galactic disk; these were in the range 
0.5-3.2xl0~ 16 s~' whenever H3 was detected (although smaller 
values were allowed but not required in cases where HJ was not 
detected). A similar estimate {(h ~ 1.2 X 10~ 16 s~') was inferred 
by Le Petit, Roueff & Herbst (2004) from observations of H+ 
and other species toward ( Persei. 
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Fig. B.l. Reduced x 1 f° r the best fit, as a function of the 
frequency assumed for the strongest H_20 + hyperfine compo- 
nent. Dotted lines indicate the values given by Miirtz et al. 
1998 and Ossenkopf et al. (2010) (P=predicted value based 
upon Strahan et al. 1986 spectroscopy; and A=astronomicalry- 
determined value) 



Appendix B: H 2 + spectroscopy 

The rest frequencies for transitions of H_20 + have been the sub- 
ject of some controversy (Ossenkopf et al. 2010). The frequency 
given by Miirtz et al. (1998) for the strongest hyperfine com- 
ponent is 1115.204 + 0.002 GHz. This value, which is based 
upon laser magnetic resonance (LMR) experiments, is ~ 28 MHz 
higher than that predicted by Ossenkopf et al. (2010) from an 
analysis of earlier LMR data (Strahan et al. 1986), and 42 MHz 
higher than the "astronomically-determined" rest frequency pre- 
sented by Ossenkopf et al. We have investigated how well the 
OH + and HaO + spectra obtained toward W49N can be recon- 
ciled for various assumed values of the H20 + transition fre- 
quency. We allowed the latter to vary, in performing the spec- 
tral fitting procedure described in §3 above, but without chang- 
ing the assumed spacing of the hyperfine components. Our best 
fit frequency for the strongest H20 + hyperfine component was 
1115.209 GHz, a value lying close to that given by Miirtz et 
al. (1998). This result would be in error if the OH + molecules 
have systematic velocity shift relative to H20 + within a given 
cloud. Also, insofar as the frequency scale is tied to the assumed 
OH + rest frequency, the fractional error in the derived H20 + fre- 
quency can be no smaller than that for the OH + rest frequency. 
Thus, our "astrophysical determination" of the H20 + rest fre- 
quency is entirely consistent with the laboratory value reported 
by Miirtz et al. (1998). However, it is apparently inconsistent 
with the significantly smaller values given by Ossenkopf et al. 
(2010). In Figure B.l, we show the reduced^ 2 for the best fit 
that can be obtained for a given assumed H20 + rest frequency, 
as a function of that frequency. The best fit shows a minimum 
close to the Miirtz et al. (1998) rest frequency and is significantly 
worse for the other values. 
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